\(~\)
\(~\)
library(Seurat)
so = Read10X(data.dir = "/mnt/market/anclab-rstudio-server/home/swil0005/projects/Wnt11/raw_data/cellranger_mRNAcounts_asgz/")
so = CreateSeuratObject(counts = so , project = "Wnt11KO" , min.cells = 20)
so
## An object of class Seurat
## 17039 features across 20586 samples within 1 assay
## Active assay: RNA (17039 features, 0 variable features)
\(~\)
\(~\)
so[["percent.mt"]] <- PercentageFeatureSet(so, pattern = "^mt-")
quantile(so@meta.data$nCount_RNA , probs = seq(0, 1, by= 0.1))
## 0% 10% 20% 30% 40% 50% 60% 70%
## 1047.0 1628.5 2651.0 3811.0 5708.0 7969.0 11349.0 14558.0
## 80% 90% 100%
## 18653.0 26267.5 200631.0
quantile(so@meta.data$nFeature_RNA , probs = seq(0, 1, by= 0.1))
## 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100%
## 242.0 631.0 704.0 889.0 1417.0 2326.0 3190.0 3670.0 4180.0 4954.5 9697.0
quantile(so@meta.data$percent.mt , probs = seq(0, 1, by= 0.1))
## 0% 10% 20% 30% 40% 50%
## 0.04941627 0.40198602 0.81072654 1.61021127 1.95670073 2.19657747
## 60% 70% 80% 90% 100%
## 2.43974133 2.79772081 4.05519233 12.10581643 38.78822197
so <- subset(so, subset = nCount_RNA > 1500 & nCount_RNA < 20000 & nFeature_RNA > 500 & nFeature_RNA < 5000 & percent.mt < 5)
so[["RNA"]]
## Assay data with 17039 features for 11686 cells
## First 10 features:
## Xkr4, Sox17, Mrpl15, Lypla1, Tcea1, Atp6v1h, Rb1cc1, 4732440D04Rik,
## Alkal1, Pcmtd1
\(~\)
\(~\)
so <- NormalizeData(so , normalization.method = "LogNormalize" , scale.factor = 10000)
\(~\)
\(~\)
so <- FindVariableFeatures(so, selection.method = "vst", nfeatures = 2000)
LabelPoints(plot = VariableFeaturePlot(so), points = head(VariableFeatures(so), 10), repel = TRUE)
\(~\)
\(~\)
library(biomaRt)
s.genes = cc.genes$s.genes
g2m.genes = cc.genes$g2m.genes
# https://rjbioinformatics.com/2016/10/14/converting-mouse-to-human-gene-names-with-biomart-package/
human = useMart("ensembl", dataset = "hsapiens_gene_ensembl")
mouse = useMart("ensembl", dataset = "mmusculus_gene_ensembl")
s.genes = getLDS(attributes = c("hgnc_symbol"), filters = "hgnc_symbol",
values = s.genes , mart = human,
attributesL = c("mgi_symbol"),martL = mouse, uniqueRows=T)
s.genes = s.genes$MGI.symbol
g2m.genes = getLDS(attributes = c("hgnc_symbol"), filters = "hgnc_symbol",
values = g2m.genes , mart = human,
attributesL = c("mgi_symbol"),martL = mouse, uniqueRows=T)
g2m.genes = g2m.genes$MGI.symbol
# mgi_symbol: gene symboles for mouse
# hgnc_symbol: gene symbols for human
\(~\)
\(~\)
# https://satijalab.org/seurat/archive/v3.2/cell_cycle_vignette.html
so <- ScaleData(so, features = rownames(so))
so = CellCycleScoring(so, s.features = s.genes, g2m.features = g2m.genes, set.ident = TRUE)
head(so[[]])
## orig.ident nCount_RNA nFeature_RNA percent.mt S.Score
## AAACCCAAGGTCGTGA-1 Wnt11KO 7798 726 0.3334188 -0.02714792
## AAACCCACAACAAGAT-1 Wnt11KO 11151 3333 2.4123397 -0.12420683
## AAACCCACAGCTCGGT-1 Wnt11KO 1572 881 1.9083969 0.12395335
## AAACCCACATGAAGCG-1 Wnt11KO 4998 683 0.4201681 -0.04567019
## AAACCCAGTGTCCGTG-1 Wnt11KO 16636 3618 2.1399375 0.05585819
## AAACCCAGTGTGTACT-1 Wnt11KO 5014 862 0.3191065 0.04130828
## G2M.Score Phase old.ident
## AAACCCAAGGTCGTGA-1 0.003610455 G2M Wnt11KO
## AAACCCACAACAAGAT-1 -0.244799795 G1 Wnt11KO
## AAACCCACAGCTCGGT-1 0.280071884 G2M Wnt11KO
## AAACCCACATGAAGCG-1 -0.013205415 G1 Wnt11KO
## AAACCCAGTGTCCGTG-1 0.484910229 G2M Wnt11KO
## AAACCCAGTGTGTACT-1 0.020404714 S Wnt11KO
RidgePlot(so, features = c("Pcna", "Top2a", "Mcm6", "Mki67"), ncol = 2)
so <- RunPCA(so, features = VariableFeatures(so), nfeatures.print = 50)
DimPlot(so)
# Cells separate by cell-cycle phase
so <- RunPCA(so, features = c(s.genes, g2m.genes) , nfeatures.print = 10)
DimPlot(so)
# so <- ScaleData(so, vars.to.regress = c("S.Score", "G2M.Score"), features = rownames(so))
so = readRDS("Wnt11KO.rds")
# Now, a PCA on the variable genes no longer returns components associated with cell cycle
so <- RunPCA(so, features = VariableFeatures(so), nfeatures.print = 50)
DimPlot(so)
# When running a PCA on only cell cycle genes, cells no longer separate by cell-cycle phase
so <- RunPCA(so, features = c(s.genes, g2m.genes))
DimPlot(so)
\(~\)
\(~\)
\(~\)
\(~\)
so = RunPCA(object = so , assay = "RNA" , features = NULL , npcs = 20 , rev.pca = FALSE ,
verbose = TRUE, ndims.print = 1:3, nfeatures.print = 10, reduction.name = "pca",
reduction.key = "PC_", seed.use = 42, approx = T)
for(i in 1:20){
DimHeatmap(so, dims = i, nfeatures = 40 , cells = 500 , reduction = "pca" ,
disp.min = -6 , disp.max = 6 , balanced = T , projected = F, ncol = NULL,
fast = TRUE, raster = TRUE, slot = "scale.data", assays = "RNA", combine = T)
}
# Selected PCs for clustering
# PCs: 1,2,3,4,5,6,7,8,13
ElbowPlot(so, ndims = 20, reduction = "pca")
\(~\)
\(~\)
so1 <- FindNeighbors(so, return.neighbor = FALSE, compute.SNN = T,
prune.SNN = 1/15, nn.method = "annoy", n.trees = 200, annoy.metric = "euclidean",
dims = c(1,2,3,4,5,6,7,8,13) , reduction = "pca" , verbose = F , do.plot = F )
# Cluster Determination
so1 = FindClusters(so1, modularity.fxn = 1 ,
initial.membership = NULL, node.sizes = NULL,
resolution = 1.5 , algorithm = 1,
n.start = 100, n.iter = 100, random.seed = 0,
group.singletons = T , temp.file.location = NULL,
edge.file.name = NULL, verbose = F)
DimPlot(so1, reduction = "pca", label = T , dims = 1:2)
so1 = RunUMAP(so1, dims = 1:20)
head(so1[[]])
## orig.ident nCount_RNA nFeature_RNA percent.mt S.Score
## AAACCCAAGGTCGTGA-1 Wnt11KO 7798 726 0.3334188 -0.03079170
## AAACCCACAACAAGAT-1 Wnt11KO 11151 3333 2.4123397 -0.11975109
## AAACCCACAGCTCGGT-1 Wnt11KO 1572 881 1.9083969 0.11639918
## AAACCCACATGAAGCG-1 Wnt11KO 4998 683 0.4201681 -0.04490011
## AAACCCAGTGTCCGTG-1 Wnt11KO 16636 3618 2.1399375 0.04813477
## AAACCCAGTGTGTACT-1 Wnt11KO 5014 862 0.3191065 0.03657714
## G2M.Score Phase old.ident RNA_snn_res.1.5 seurat_clusters
## AAACCCAAGGTCGTGA-1 -0.008638847 G1 Wnt11KO 2 2
## AAACCCACAACAAGAT-1 -0.257422641 G1 Wnt11KO 7 7
## AAACCCACAGCTCGGT-1 0.263503318 G2M Wnt11KO 8 8
## AAACCCACATGAAGCG-1 -0.040737166 G1 Wnt11KO 1 1
## AAACCCAGTGTCCGTG-1 0.472502356 G2M Wnt11KO 14 14
## AAACCCAGTGTGTACT-1 0.005532894 S Wnt11KO 1 1
names(so1)
## [1] "RNA" "RNA_nn" "RNA_snn" "pca" "umap"
DimPlot(so1, reduction = "umap", label = T)
\(~\)
\(~\)
features = "Six2"
VlnPlot(so1, features = features) + NoLegend()
RidgePlot(object = so1 , features = features) + NoLegend()
DotPlot(so1, features = features) + RotatedAxis()
FeaturePlot(so1, features = features)
\(~\)
\(~\)
\(~\)
Wnt4, C1qtnf12 (Fam132a), and Pax8 are the markers for Committing NP (N_Com).
S100g, Mal, Wfdc2, Aqp3, and Aqp2 are the markers for Distal nephron (N_Dist).
Wnt9b and Aqp2 are the markers for UE cortical collecting duct (UE_Ccd).
Upk3a, Upk1a, and Upk1b are the markers for UE Medullary collecting duct and Urothelium (UE_Mcd-Uro).
features = "Wnt4"
VlnPlot(so1, features = features) + NoLegend()
DotPlot(so1, features = features) + RotatedAxis()
FeaturePlot(so1, features = features)
features = "S100g"
VlnPlot(so1, features = features) + NoLegend()
DotPlot(so1, features = features) + RotatedAxis()
FeaturePlot(so1, features = features)
features = "Wnt9b"
VlnPlot(so1, features = features) + NoLegend()
DotPlot(so1, features = features) + RotatedAxis()
FeaturePlot(so1, features = features)
features = "Upk3a"
VlnPlot(so1, features = features) + NoLegend()
DotPlot(so1, features = features) + RotatedAxis()
FeaturePlot(so1, features = features)
\(~\)
\(~\)
### Nephron
# Nephron progenitor (N_NP), (Cited1, Crym, Six2 in clusters 4, 19)
# Committing NP/Committing nephron (PTA/RV/SSB)(N_Com), (Wnt4, in clusters 17,24), (C1qtnf12, inclusters 19,24,25,26), (Lef1 in cluster 17), (Pax8 in 24,26)
# Distal nephron/connecting segment (N_Dist), (S100g in 13,24,26), (Mal,Wfdc2,Aqp3,Aqp2 in 13,24), (Slc12a1 in 13)
# Nephron medial/Loop of Henle (N_LOH), No cluster.
# marker of Nephron proximal (N_Prox), (Hnf4a, Lrp2, Aldob in 26)
# Nephron podocyte (N_Pod), (Nphs1, Nphs2, Mafb, Synpo in 18)
### Ureteric Epithelium (Calb1 is the marker upregulated in 13)
# UE tip (UE_Tip), (Ret,Wnt11,Crlf1 in 13).
# UE cortical collecting duct (UE_Ccd), (Wnt9b,Aqp2 in 13, 24), (Slit2 in 13).
# UE Medullary collecting duct and Urothelium (UE_Mcd-Uro), (Upk3a, Upk1a in 24), (Upk1b in 24,25).
### Stroma (Pdgfra, Pdgfrb in 3,6,7,9,11,12,15,17,21)
# Cortical stroma (S_CS), (Ntn1 in 6,12,21), (Foxd1 in 6,11,12,15,21),
# Collecting-duct associated (S_CD), (Lef1 in 9,17), (Alx1 in 9,17)
# Stroma surrounding ureter S_UR, (Tbx18 in 11),
# Generic stromal markers S,
### Endothelium
# Endothelial/vascular (Endo), (Pecam1 in clusters 20,23), (Cdh5, Flt1, Kdr in 20).
### Immune
# Immune cells/myleoid cells (Imm), (Lyz2 in 10,14 (main),16,22,23),
# (C1qa in 10,14(a little),16,22), (Fcer1g in 10, 14, 16, 22, 23)
\(~\)
\(~\)
Current.cluster.ids <- c(0:26)
New.cluster.ids <- c("Blood", "Blood", "Blood", "Stroma", "N_NP", "Blood", "S_CS" , "Stroma" , "8" ,
"S_CD" , "Imm" , "S_UR" , "S_CS" , "N_Dist/UE_Tip/UE_Ccd" , "Imm" , "S_CS" ,
"Imm" , "S_CD" , "N_Pod" , "N_NP" , "Endo" , "S_CS" , "Imm" , "Imm/Endo" ,
"N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd" , "N_Com/UE_Mcd-Uro" , "N_Prox/N_Com/N_Dist")
names(New.cluster.ids) <- levels(so1)
so1 <- RenameIdents(so1, New.cluster.ids)
DimPlot(so1, reduction = "umap", label = TRUE, pt.size = 0.5 , label.size = 3) + NoLegend()
\(~\)
so1 = AddMetaData(so1 , Idents(so1) , "Cell_Type")
head(so1[[]] , 20)
## orig.ident nCount_RNA nFeature_RNA percent.mt S.Score
## AAACCCAAGGTCGTGA-1 Wnt11KO 7798 726 0.3334188 -0.03079170
## AAACCCACAACAAGAT-1 Wnt11KO 11151 3333 2.4123397 -0.11975109
## AAACCCACAGCTCGGT-1 Wnt11KO 1572 881 1.9083969 0.11639918
## AAACCCACATGAAGCG-1 Wnt11KO 4998 683 0.4201681 -0.04490011
## AAACCCAGTGTCCGTG-1 Wnt11KO 16636 3618 2.1399375 0.04813477
## AAACCCAGTGTGTACT-1 Wnt11KO 5014 862 0.3191065 0.03657714
## AAACCCATCAACACGT-1 Wnt11KO 15992 4076 3.0077539 -0.20074218
## AAACCCATCCGCAAAT-1 Wnt11KO 7300 724 0.2739726 0.01221797
## AAACCCATCTATCGTT-1 Wnt11KO 7010 713 1.0271041 0.03923562
## AAACCCATCTTGGCTC-1 Wnt11KO 1696 1023 2.2405660 0.07165411
## AAACGAAAGACCTGGA-1 Wnt11KO 12338 3359 1.9289998 -0.07249333
## AAACGAAAGCAGATAT-1 Wnt11KO 1721 806 1.1621150 -0.03952054
## AAACGAAAGCCTGACC-1 Wnt11KO 1695 744 1.2979351 0.19732527
## AAACGAAAGGAACGTC-1 Wnt11KO 2101 689 0.4759638 0.15713407
## AAACGAAAGGGCGAAG-1 Wnt11KO 7994 2793 2.4018014 -0.15401931
## AAACGAAAGGTCGTGA-1 Wnt11KO 12729 2934 2.3803912 -0.03005646
## AAACGAAAGTCATGAA-1 Wnt11KO 16758 3936 2.1780642 -0.13775837
## AAACGAACAAACTAAG-1 Wnt11KO 15533 3460 1.5386596 0.25344568
## AAACGAACAATGTCTG-1 Wnt11KO 3199 1473 1.1878712 0.20664185
## AAACGAAGTCCAGAAG-1 Wnt11KO 17039 3391 2.0775867 -0.19162184
## G2M.Score Phase old.ident RNA_snn_res.1.5 seurat_clusters
## AAACCCAAGGTCGTGA-1 -0.008638847 G1 Wnt11KO 2 2
## AAACCCACAACAAGAT-1 -0.257422641 G1 Wnt11KO 7 7
## AAACCCACAGCTCGGT-1 0.263503318 G2M Wnt11KO 8 8
## AAACCCACATGAAGCG-1 -0.040737166 G1 Wnt11KO 1 1
## AAACCCAGTGTCCGTG-1 0.472502356 G2M Wnt11KO 14 14
## AAACCCAGTGTGTACT-1 0.005532894 S Wnt11KO 1 1
## AAACCCATCAACACGT-1 -0.281983255 G1 Wnt11KO 17 17
## AAACCCATCCGCAAAT-1 -0.096170188 S Wnt11KO 2 2
## AAACCCATCTATCGTT-1 -0.046707833 S Wnt11KO 0 0
## AAACCCATCTTGGCTC-1 -0.244101503 S Wnt11KO 8 8
## AAACGAAAGACCTGGA-1 0.139565277 G2M Wnt11KO 12 12
## AAACGAAAGCAGATAT-1 0.016302728 G2M Wnt11KO 1 1
## AAACGAAAGCCTGACC-1 0.038377154 S Wnt11KO 8 8
## AAACGAAAGGAACGTC-1 -0.010733195 S Wnt11KO 1 1
## AAACGAAAGGGCGAAG-1 -0.313363478 G1 Wnt11KO 21 21
## AAACGAAAGGTCGTGA-1 0.432119583 G2M Wnt11KO 20 20
## AAACGAAAGTCATGAA-1 -0.194367466 G1 Wnt11KO 13 13
## AAACGAACAAACTAAG-1 -0.191072581 S Wnt11KO 16 16
## AAACGAACAATGTCTG-1 -0.074074612 S Wnt11KO 8 8
## AAACGAAGTCCAGAAG-1 -0.295034457 G1 Wnt11KO 8 8
## Cell_Type
## AAACCCAAGGTCGTGA-1 Blood
## AAACCCACAACAAGAT-1 Stroma
## AAACCCACAGCTCGGT-1 8
## AAACCCACATGAAGCG-1 Blood
## AAACCCAGTGTCCGTG-1 Imm
## AAACCCAGTGTGTACT-1 Blood
## AAACCCATCAACACGT-1 S_CD
## AAACCCATCCGCAAAT-1 Blood
## AAACCCATCTATCGTT-1 Blood
## AAACCCATCTTGGCTC-1 8
## AAACGAAAGACCTGGA-1 S_CS
## AAACGAAAGCAGATAT-1 Blood
## AAACGAAAGCCTGACC-1 8
## AAACGAAAGGAACGTC-1 Blood
## AAACGAAAGGGCGAAG-1 S_CS
## AAACGAAAGGTCGTGA-1 Endo
## AAACGAAAGTCATGAA-1 N_Dist/UE_Tip/UE_Ccd
## AAACGAACAAACTAAG-1 Imm
## AAACGAACAATGTCTG-1 8
## AAACGAAGTCCAGAAG-1 8
\(~\)
\(~\)
so.markers <- FindAllMarkers(so1, slot = "data", min.pct = 0.1, verbose = F ,
logfc.threshold = 0.5 , test.use = "t", assay = "RNA",
min.cells.feature = 10 , min.cells.group = 10)
# All Markers with absolute average log fold change larger than two
subset(so.markers , abs(so.markers$avg_log2FC) > 2 & so.markers$p_val_adj < 0.05)[,c(2,5,6)]
## avg_log2FC p_val_adj cluster
## Hba-a1 5.009865 0.000000e+00 Blood
## Hba-a2 5.001911 0.000000e+00 Blood
## Hbb-bt 4.938371 0.000000e+00 Blood
## Hbb-bs 4.919223 0.000000e+00 Blood
## Alas2 3.486368 0.000000e+00 Blood
## Mkrn1 2.993951 0.000000e+00 Blood
## Bpgm 2.690968 0.000000e+00 Blood
## Snca 2.677801 0.000000e+00 Blood
## Fech 2.190977 0.000000e+00 Blood
## Fam220a 2.070838 0.000000e+00 Blood
## Slc25a37 2.019208 0.000000e+00 Blood
## Cfl1 -2.042088 0.000000e+00 Blood
## Col3a1 -2.048892 0.000000e+00 Blood
## Lgals1 -2.090133 0.000000e+00 Blood
## Col1a2 -2.109223 0.000000e+00 Blood
## Meg3 -2.155789 0.000000e+00 Blood
## Vim -2.172572 0.000000e+00 Blood
## Sparc -2.204993 0.000000e+00 Blood
## Cst3 -2.211365 0.000000e+00 Blood
## Itm2b -2.216454 0.000000e+00 Blood
## Actb -2.217154 0.000000e+00 Blood
## Dcn -2.246423 0.000000e+00 Blood
## Malat1 -2.263107 0.000000e+00 Blood
## mt-Cytb -2.264501 0.000000e+00 Blood
## Actg1 -2.291778 0.000000e+00 Blood
## Col1a1 -2.350176 0.000000e+00 Blood
## Cdkn1c -2.438399 0.000000e+00 Blood
## Tmsb4x -2.997531 0.000000e+00 Blood
## Apoe -3.332003 4.147493e-249 Blood
## Fcer1g -2.276320 2.509331e-195 Blood
## Tyrobp -2.213460 8.715077e-191 Blood
## Hbb-y 4.038364 1.749192e-137 Blood
## Lyz2 -2.411889 1.689527e-133 Blood
## C1qb -2.222077 6.594791e-130 Blood
## Hba-x 3.539685 4.496074e-62 Blood
## Mkrn1.1 -2.232286 0.000000e+00 Stroma
## Alas2.1 -2.541719 0.000000e+00 Stroma
## Hbb-bt.1 -4.214194 0.000000e+00 Stroma
## Hbb-bs.1 -4.319514 0.000000e+00 Stroma
## Hba-a2.1 -4.396159 0.000000e+00 Stroma
## Hba-a1.1 -4.424480 0.000000e+00 Stroma
## Fcer1g.1 -2.003609 2.477138e-166 Stroma
## Lyz2.1 -2.176337 2.058066e-146 Stroma
## Mgp.1 2.009866 2.396092e-121 Stroma
## Hbb-y.1 -4.876611 9.283181e-85 Stroma
## Apoe.1 -2.749975 9.319485e-55 Stroma
## Hba-x.1 -3.763994 9.329004e-50 Stroma
## Fxyd6.1 2.077086 0.000000e+00 N_NP
## Col1a2.2 -2.016121 0.000000e+00 N_NP
## Col3a1.2 -2.087078 0.000000e+00 N_NP
## Col1a1.2 -2.210267 0.000000e+00 N_NP
## Alas2.2 -2.618031 0.000000e+00 N_NP
## Tmsb4x.1 -3.172623 0.000000e+00 N_NP
## Hbb-bt.2 -4.557309 0.000000e+00 N_NP
## Hbb-bs.2 -4.608028 0.000000e+00 N_NP
## Hba-a2.2 -4.699828 0.000000e+00 N_NP
## Hba-a1.2 -4.783347 0.000000e+00 N_NP
## Cited1.2 2.723417 1.683677e-319 N_NP
## Mkrn1.2 -2.099821 4.108987e-279 N_NP
## Apoe.2 -3.369404 3.731432e-243 N_NP
## Fcer1g.2 -2.043260 1.830185e-202 N_NP
## Dcn.2 -2.064448 1.588360e-199 N_NP
## Lyz2.2 -2.176970 3.892392e-171 N_NP
## C1qb.2 -2.018274 1.843467e-153 N_NP
## Hbb-y.2 -4.665824 4.162694e-94 N_NP
## Hba-x.2 -3.730137 5.668646e-41 N_NP
## Foxd1.1 2.210584 0.000000e+00 S_CS
## Aldh1a2.2 2.079697 0.000000e+00 S_CS
## Mkrn1.3 -2.214426 0.000000e+00 S_CS
## Alas2.3 -2.540206 0.000000e+00 S_CS
## Hbb-bs.3 -4.185937 0.000000e+00 S_CS
## Hbb-bt.3 -4.197440 0.000000e+00 S_CS
## Hba-a2.3 -4.351547 0.000000e+00 S_CS
## Hba-a1.3 -4.356825 0.000000e+00 S_CS
## Meg3.2 2.070969 2.825017e-251 S_CS
## Dlk1.2 2.285813 1.035762e-245 S_CS
## Tyrobp.3 -2.030311 2.367143e-226 S_CS
## Fcer1g.3 -2.078628 8.674893e-217 S_CS
## Lyz2.3 -2.192568 3.745981e-157 S_CS
## C1qb.3 -2.013512 1.427166e-135 S_CS
## Hbb-y.3 -3.257407 1.847239e-53 S_CS
## Apoe.3 -2.795011 1.748306e-49 S_CS
## Hba-x.3 -2.853523 2.040742e-23 S_CS
## Cdkn1c.2 2.686545 4.726188e-293 S_CD
## Alas2.5 -2.389891 1.377525e-245 S_CD
## Hba-a1.4 -4.034875 5.933338e-218 S_CD
## Hba-a2.4 -4.042353 1.083976e-217 S_CD
## Hbb-bt.5 -3.904486 2.539240e-217 S_CD
## Hbb-bs.5 -3.833715 2.343454e-216 S_CD
## Penk.3 2.604979 8.623620e-178 S_CD
## Fcer1g.5 -2.023176 6.655842e-177 S_CD
## Lyz2.5 -2.120207 9.388315e-114 S_CD
## Hbb-y.4 -3.490904 1.661583e-27 S_CD
## Apoe.5 -2.685152 4.851766e-23 S_CD
## Hba-x.4 -2.755749 4.275659e-11 S_CD
## Apoe.6 5.745254 0.000000e+00 Imm
## Lyz2.6 4.898655 0.000000e+00 Imm
## C1qb.6 4.811275 0.000000e+00 Imm
## Fcer1g.6 4.660813 0.000000e+00 Imm
## Tyrobp.6 4.386300 0.000000e+00 Imm
## Coro1a.6 3.587922 0.000000e+00 Imm
## Cd52.6 3.573172 0.000000e+00 Imm
## Ctss.6 3.500510 0.000000e+00 Imm
## Ctsc.6 3.261860 0.000000e+00 Imm
## Ctsb.6 3.247724 0.000000e+00 Imm
## Tmsb4x.5 3.221111 0.000000e+00 Imm
## Lst1.6 3.209396 0.000000e+00 Imm
## Ms4a6c.6 3.206885 0.000000e+00 Imm
## Cyba.3 3.084055 0.000000e+00 Imm
## Selenop.3 3.039459 0.000000e+00 Imm
## Arhgdib.6 3.018058 0.000000e+00 Imm
## Aif1.6 2.919618 0.000000e+00 Imm
## Psap.2 2.855854 0.000000e+00 Imm
## Lgmn.5 2.823587 0.000000e+00 Imm
## Ptpn18.5 2.812157 0.000000e+00 Imm
## Fcgr3.6 2.734368 0.000000e+00 Imm
## Spi1.5 2.692691 0.000000e+00 Imm
## Ctsd.1 2.684440 0.000000e+00 Imm
## Cst3.4 2.678812 0.000000e+00 Imm
## Rac2.6 2.671791 0.000000e+00 Imm
## Laptm5.5 2.609907 0.000000e+00 Imm
## Ctsh.5 2.572160 0.000000e+00 Imm
## Ly86.5 2.557617 0.000000e+00 Imm
## B2m.3 2.541774 0.000000e+00 Imm
## Sat1.5 2.537128 0.000000e+00 Imm
## Csf1r.5 2.531573 0.000000e+00 Imm
## Arpc1b.2 2.499485 0.000000e+00 Imm
## Ms4a6b.5 2.477358 0.000000e+00 Imm
## Itm2b.3 2.420265 0.000000e+00 Imm
## Fyb.5 2.390871 0.000000e+00 Imm
## Gmfg.5 2.345757 0.000000e+00 Imm
## Cd68.2 2.331817 0.000000e+00 Imm
## Ctsz.2 2.329625 0.000000e+00 Imm
## Unc93b1.2 2.302544 0.000000e+00 Imm
## H2-D1.5 2.289112 0.000000e+00 Imm
## Ckb.3 2.282515 0.000000e+00 Imm
## Grn.1 2.208881 0.000000e+00 Imm
## Hexa.1 2.175718 0.000000e+00 Imm
## Rgs10.2 2.164107 0.000000e+00 Imm
## Pld4.1 2.160643 0.000000e+00 Imm
## Sh3bgrl3.2 2.152989 0.000000e+00 Imm
## Lamp1.1 2.125554 0.000000e+00 Imm
## Cotl1.2 2.108484 0.000000e+00 Imm
## Clta.1 2.095699 0.000000e+00 Imm
## Clec4a2.1 2.081758 0.000000e+00 Imm
## Ehd4.1 2.078501 0.000000e+00 Imm
## Tspo.3 2.061047 0.000000e+00 Imm
## Psmb8.1 2.034493 0.000000e+00 Imm
## Mdk.4 -2.186830 0.000000e+00 Imm
## Sparc.3 -2.251869 0.000000e+00 Imm
## Cdkn1c.3 -2.697607 0.000000e+00 Imm
## C1qc.6 4.270142 1.168472e-316 Imm
## Col1a1.5 -2.105762 1.330768e-311 Imm
## Cx3cr1.1 2.239852 2.777859e-306 Imm
## Hba-a1.5 -3.602695 5.153050e-292 Imm
## Hba-a2.5 -3.576468 3.593782e-286 Imm
## Hbb-bs.6 -3.503388 3.790010e-284 Imm
## Alas2.6 -2.255738 1.697110e-282 Imm
## C1qa.6 4.337589 3.379972e-280 Imm
## Hbb-bt.6 -3.384944 2.564061e-268 Imm
## Alox5ap.6 2.438275 1.453811e-256 Imm
## Slfn2.1 2.004201 6.142255e-255 Imm
## Dcn.4 -2.023512 1.164122e-211 Imm
## Ly6e.2 2.001744 5.593573e-209 Imm
## Fos.4 2.541716 2.365344e-203 Imm
## Fcrls.1 2.295141 8.592411e-190 Imm
## Mt1.5 2.158338 4.360272e-183 Imm
## Dab2.2 2.155805 2.597052e-180 Imm
## Pf4.6 4.001026 2.977531e-172 Imm
## Ifitm3.3 2.782791 1.549683e-171 Imm
## Ms4a7.1 2.081041 3.561917e-163 Imm
## Mrc1 2.224318 7.683972e-145 Imm
## Lgals3.1 2.463998 1.186633e-125 Imm
## Ifi27l2a.1 2.220170 3.375187e-124 Imm
## Hmox1.3 2.457087 1.475244e-106 Imm
## Cd74 2.227840 3.547063e-52 Imm
## Hbb-y.5 -4.356507 7.748176e-48 Imm
## Hba-x.5 -3.434211 2.457970e-32 Imm
## Col1a1.6 2.963275 0.000000e+00 S_UR
## Col1a2.6 2.604663 0.000000e+00 S_UR
## Sparc.4 2.369315 0.000000e+00 S_UR
## Col3a1.6 2.552808 2.900421e-279 S_UR
## Dcn.5 3.195732 2.058517e-259 S_UR
## Meg3.5 2.254437 7.982242e-175 S_UR
## Alas2.7 -2.412515 4.183164e-168 S_UR
## Akap12.6 2.245148 1.563906e-158 S_UR
## Cxcl12.4 2.371043 5.286460e-147 S_UR
## Lum.4 3.188742 1.066729e-145 S_UR
## Hbb-bs.7 -3.954088 3.901525e-144 S_UR
## Hba-a1.6 -4.069833 6.856096e-142 S_UR
## Hbb-bt.7 -3.921396 2.382666e-140 S_UR
## Dlk1.6 2.247238 8.160804e-140 S_UR
## Hba-a2.6 -4.050937 1.924021e-139 S_UR
## Sfrp2.3 2.312653 6.549013e-134 S_UR
## Mkrn1.7 -2.098201 4.846172e-127 S_UR
## Lyz2.7 -2.075025 1.825218e-77 S_UR
## Hbb-y.6 -2.502790 6.691056e-12 S_UR
## Hba-x.6 -2.829514 5.347443e-06 S_UR
## Epcam.7 3.468537 3.946723e-312 N_Dist/UE_Tip/UE_Ccd
## Wfdc2.6 3.437953 1.338285e-270 N_Dist/UE_Tip/UE_Ccd
## Ldhb.2 2.238744 1.157932e-254 N_Dist/UE_Tip/UE_Ccd
## Cldn7.1 2.586154 5.531850e-247 N_Dist/UE_Tip/UE_Ccd
## Vim.6 -2.410637 2.872840e-222 N_Dist/UE_Tip/UE_Ccd
## Col1a1.7 -2.269590 7.278929e-183 N_Dist/UE_Tip/UE_Ccd
## Col3a1.7 -2.117079 1.399181e-182 N_Dist/UE_Tip/UE_Ccd
## Col1a2.7 -2.050271 8.968457e-175 N_Dist/UE_Tip/UE_Ccd
## Krt18.6 2.524894 1.911047e-173 N_Dist/UE_Tip/UE_Ccd
## Lgals1.6 -2.364860 1.860429e-168 N_Dist/UE_Tip/UE_Ccd
## Calb1.7 4.093941 7.014694e-165 N_Dist/UE_Tip/UE_Ccd
## Dcn.6 -2.162377 2.647358e-159 N_Dist/UE_Tip/UE_Ccd
## Alas2.8 -2.471052 1.317094e-158 N_Dist/UE_Tip/UE_Ccd
## Cldn3 2.374805 8.619318e-156 N_Dist/UE_Tip/UE_Ccd
## Krt8.3 2.192433 7.328108e-155 N_Dist/UE_Tip/UE_Ccd
## Hbb-bt.8 -4.551041 9.181954e-149 N_Dist/UE_Tip/UE_Ccd
## Hba-a1.7 -4.492709 1.992368e-147 N_Dist/UE_Tip/UE_Ccd
## Hbb-bs.8 -4.219864 3.184453e-147 N_Dist/UE_Tip/UE_Ccd
## Hba-a2.7 -4.490026 2.130870e-146 N_Dist/UE_Tip/UE_Ccd
## Rprm 2.748766 1.448208e-124 N_Dist/UE_Tip/UE_Ccd
## Cdkn1c.5 -2.297807 1.910830e-98 N_Dist/UE_Tip/UE_Ccd
## Lyz2.8 -2.066784 1.228262e-62 N_Dist/UE_Tip/UE_Ccd
## Apoe.8 -2.883635 4.703977e-45 N_Dist/UE_Tip/UE_Ccd
## Hbb-y.7 -3.422496 2.162567e-33 N_Dist/UE_Tip/UE_Ccd
## Hba-x.7 -2.642879 3.408679e-09 N_Dist/UE_Tip/UE_Ccd
## S100g 2.047456 2.851832e-08 N_Dist/UE_Tip/UE_Ccd
## Igfbp7.5 2.515136 5.162524e-11 N_Pod
## Hba-a2.8 -2.107410 2.210355e-10 N_Pod
## Hba-a1.8 -2.074964 3.090650e-10 N_Pod
## Mafb.2 3.063028 3.590388e-10 N_Pod
## Hbb-bs.9 -2.002717 1.159542e-09 N_Pod
## Pcp4 2.024515 1.121237e-08 N_Pod
## Rasl11a.1 2.138632 2.426344e-08 N_Pod
## R3hdml 2.656770 3.874996e-06 N_Pod
## Apoe.9 -2.388812 4.481074e-02 N_Pod
## Egfl7 3.051555 8.751735e-65 Endo
## Crip2.2 2.349887 3.234954e-56 Endo
## Fabp5.2 2.753293 4.083615e-56 Endo
## Fabp4 4.609578 7.462791e-54 Endo
## Icam2 2.397084 1.298368e-52 Endo
## Emcn 2.894633 3.426105e-52 Endo
## Cldn5.1 2.632344 4.089893e-49 Endo
## Plvap 3.200724 1.884233e-47 Endo
## Ramp2.4 2.269281 6.870409e-47 Endo
## Gng11.2 2.163268 1.007255e-43 Endo
## Fam167b 2.115451 3.655401e-43 Endo
## Lyz2.10 -2.007078 3.782979e-18 Endo
## Hba-a1.9 -2.301117 3.889220e-18 Endo
## Hba-a2.9 -2.273174 8.589599e-18 Endo
## Hbb-bt.10 -2.367982 1.790171e-17 Endo
## Hbb-bs.10 -2.128467 1.834480e-17 Endo
## Apoe.10 -2.178600 2.360740e-04 Endo
## Coro1a.11 2.731706 8.312548e-61 Imm/Endo
## Arhgdib.9 2.809156 7.675932e-59 Imm/Endo
## Rac2.11 3.068671 2.047729e-53 Imm/Endo
## Sh3bgrl3.5 2.053310 1.830868e-49 Imm/Endo
## Gmfg.9 2.396527 1.660922e-37 Imm/Endo
## Srgn.4 3.996566 1.586487e-35 Imm/Endo
## Tyrobp.11 3.063740 1.431506e-34 Imm/Endo
## Fcer1g.11 2.031083 2.911497e-34 Imm/Endo
## Laptm5.9 2.037861 2.497048e-31 Imm/Endo
## Tspo.5 2.046180 7.119600e-31 Imm/Endo
## Cd52.11 2.795578 2.078738e-30 Imm/Endo
## H2-D1.8 2.372063 3.636345e-30 Imm/Endo
## Lcp1.2 2.217571 1.629972e-28 Imm/Endo
## Alas2.11 -2.261315 3.901658e-25 Imm/Endo
## Fxyd5.2 3.127217 5.593747e-25 Imm/Endo
## Man2b1.2 2.108062 1.040911e-19 Imm/Endo
## Emb.2 2.007622 1.244493e-19 Imm/Endo
## Hba-a2.10 -3.359457 1.382025e-19 Imm/Endo
## Fcgr3.11 2.035771 3.357833e-19 Imm/Endo
## Hba-a1.10 -3.370828 4.082851e-19 Imm/Endo
## Hbb-bt.11 -3.378443 1.510421e-18 Imm/Endo
## Ptprcap 2.320070 5.492539e-18 Imm/Endo
## Hbb-bs.11 -3.134245 8.107288e-18 Imm/Endo
## Alox5ap.11 3.542391 2.120266e-17 Imm/Endo
## Hdc 3.564075 4.998452e-14 Imm/Endo
## Ifitm1.5 6.049862 1.366749e-13 Imm/Endo
## Cd9.1 2.033821 2.527257e-13 Imm/Endo
## Ltb 2.184038 8.656467e-13 Imm/Endo
## Plac8.6 2.310042 1.574734e-12 Imm/Endo
## Lgals3.2 2.276612 4.565328e-12 Imm/Endo
## Pglyrp1 2.968278 9.102910e-10 Imm/Endo
## Cpa3 4.401511 2.154801e-09 Imm/Endo
## Cyp11a1 2.382940 8.700792e-09 Imm/Endo
## Cd69 2.116353 2.276980e-07 Imm/Endo
## Cd200r3 2.171432 4.762573e-07 Imm/Endo
## Ccl9.1 2.723936 1.054960e-06 Imm/Endo
## Ccl3.1 2.599118 2.081037e-05 Imm/Endo
## Mcpt8 3.376087 4.712338e-04 Imm/Endo
## Krt18.8 4.209438 3.437371e-80 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Epcam.11 3.125237 1.942720e-78 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Krt8.4 4.045806 1.024845e-73 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Alas2.12 -2.661593 5.231427e-71 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Cldn7.3 3.280594 4.502295e-69 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Spint2.2 2.404951 9.420929e-66 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Wfdc2.9 3.015752 1.573198e-60 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Tacstd2.1 3.422996 4.103121e-53 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Bcam.2 2.469606 3.258920e-51 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Col1a2.11 -2.114803 4.988994e-50 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Krt7.2 3.319589 1.572614e-48 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Atp1b1.3 2.264765 3.020460e-48 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Hbb-bt.12 -5.099266 7.968804e-47 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Lgals1.10 -2.614273 1.016680e-44 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Vim.8 -2.371101 1.127828e-44 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Hba-a1.11 -4.945765 1.588618e-44 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Krt19 3.707714 6.416738e-44 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Hba-a2.11 -4.982656 1.626493e-43 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Mest.10 -2.184883 2.313260e-43 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Hbb-bs.12 -4.666126 2.535880e-43 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Tspan8 2.936064 3.965445e-43 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Cldn4.1 2.244381 4.109752e-43 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Col1a1.11 -2.267846 1.027212e-42 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Lyz2.12 -2.131684 8.184368e-42 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Foxq1 2.647323 3.888820e-41 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Col3a1.11 -2.082213 6.245223e-41 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Hspb1.1 2.579797 2.757308e-38 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Mal.1 2.477846 5.447142e-38 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Cystm1 2.077164 2.978615e-37 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Wnt4.1 2.175289 5.372146e-34 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Pdzk1ip1.1 2.933248 1.845267e-32 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Dcn.10 -2.097509 1.388950e-31 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Apoe.11 -3.355543 3.290511e-28 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Cldn3.1 2.177516 1.069880e-26 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Fxyd4 3.719893 4.311397e-25 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Cdh16.2 2.275180 9.254549e-25 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Fxyd3 3.074420 6.887786e-24 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Id1.6 2.145059 1.767947e-23 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Sfn 2.090866 2.120981e-20 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## S100a6.10 2.074633 1.483915e-16 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Cdkn1c.9 -2.041002 3.322340e-13 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Upk3a 2.534381 4.385484e-12 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Hbb-y.10 -5.007761 7.326094e-09 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Hba-x.9 -3.880112 1.676136e-07 N_Com/N_Dist/UE_Mcd-Uro/UE_Ccd
## Sparc.8 3.346582 1.980169e-57 N_Com/UE_Mcd-Uro
## Col3a1.12 2.673101 2.064136e-54 N_Com/UE_Mcd-Uro
## Sfrp2.8 3.431200 4.023557e-51 N_Com/UE_Mcd-Uro
## Col1a1.12 2.735051 5.039510e-49 N_Com/UE_Mcd-Uro
## Dlk1.12 3.566384 8.619757e-46 N_Com/UE_Mcd-Uro
## Fbln2.2 3.366598 5.894116e-45 N_Com/UE_Mcd-Uro
## Cryab.1 2.946992 9.215705e-45 N_Com/UE_Mcd-Uro
## Krt19.1 3.330028 1.886138e-44 N_Com/UE_Mcd-Uro
## Upk3b 3.784039 4.428858e-42 N_Com/UE_Mcd-Uro
## Crip1.8 3.226252 5.677367e-40 N_Com/UE_Mcd-Uro
## Meg3.10 2.424259 6.181061e-39 N_Com/UE_Mcd-Uro
## Slc9a3r1.4 2.503056 6.862073e-38 N_Com/UE_Mcd-Uro
## Hspb1.2 3.140294 9.232364e-38 N_Com/UE_Mcd-Uro
## Lyz2.13 -2.137057 1.161133e-37 N_Com/UE_Mcd-Uro
## Col14a1.10 2.756124 1.945911e-37 N_Com/UE_Mcd-Uro
## Cavin3.4 2.155528 3.661515e-36 N_Com/UE_Mcd-Uro
## Col1a2.12 2.168213 4.391321e-36 N_Com/UE_Mcd-Uro
## Krt7.3 2.819083 1.586142e-35 N_Com/UE_Mcd-Uro
## S100a6.11 3.476857 2.269283e-34 N_Com/UE_Mcd-Uro
## Tmem108 2.200778 1.831667e-30 N_Com/UE_Mcd-Uro
## Itm2a.9 3.388282 1.741411e-29 N_Com/UE_Mcd-Uro
## Hbb-bt.13 -4.608916 3.334300e-28 N_Com/UE_Mcd-Uro
## Alas2.13 -2.393809 1.227427e-27 N_Com/UE_Mcd-Uro
## Rspo1.1 2.009231 1.479489e-27 N_Com/UE_Mcd-Uro
## Hbb-bs.13 -4.124014 4.626799e-27 N_Com/UE_Mcd-Uro
## Hba-a1.12 -4.389347 4.757441e-27 N_Com/UE_Mcd-Uro
## Hba-a2.12 -4.379492 6.346572e-27 N_Com/UE_Mcd-Uro
## Eln 2.327724 2.249836e-25 N_Com/UE_Mcd-Uro
## Rarres2.10 2.427102 6.673742e-23 N_Com/UE_Mcd-Uro
## Aebp1.3 2.045606 1.125318e-22 N_Com/UE_Mcd-Uro
## 1500015O10Rik.1 2.264760 2.178616e-22 N_Com/UE_Mcd-Uro
## Igfbp6.6 2.247473 2.708995e-12 N_Com/UE_Mcd-Uro
## Fmo2 3.033400 9.313511e-12 N_Com/UE_Mcd-Uro
## Hba-x.10 -3.960187 1.432245e-11 N_Com/UE_Mcd-Uro
## Hbb-y.11 -5.118240 5.472239e-10 N_Com/UE_Mcd-Uro
## Stmn2.2 3.272816 8.936980e-08 N_Com/UE_Mcd-Uro
## Cyba.5 2.354472 7.533231e-35 N_Prox/N_Com/N_Dist
## Fxyd2.1 5.291368 2.245865e-34 N_Prox/N_Com/N_Dist
## Ly6a.3 4.212719 5.478547e-33 N_Prox/N_Com/N_Dist
## Lrpap1.2 2.651403 1.108613e-31 N_Prox/N_Com/N_Dist
## Calml4.2 3.179583 9.982252e-31 N_Prox/N_Com/N_Dist
## Slc9a3r1.5 2.215087 1.068975e-30 N_Prox/N_Com/N_Dist
## Mpc2.4 2.486508 1.197684e-30 N_Prox/N_Com/N_Dist
## Pdzk1 2.746115 1.829981e-29 N_Prox/N_Com/N_Dist
## Rida.1 3.037432 1.697521e-28 N_Prox/N_Com/N_Dist
## Chchd10.4 2.554388 1.904516e-27 N_Prox/N_Com/N_Dist
## Amn 2.654958 1.551127e-26 N_Prox/N_Com/N_Dist
## Mdh1.4 2.714683 2.462354e-26 N_Prox/N_Com/N_Dist
## Igfbp7.10 2.365169 3.596736e-26 N_Prox/N_Com/N_Dist
## Folr1.2 2.530971 9.119426e-26 N_Prox/N_Com/N_Dist
## Ezr.4 2.021599 9.895440e-25 N_Prox/N_Com/N_Dist
## Atp1b1.5 2.674680 1.059714e-24 N_Prox/N_Com/N_Dist
## Akr7a5.4 2.170087 1.805664e-24 N_Prox/N_Com/N_Dist
## Npnt.4 2.169077 1.910573e-24 N_Prox/N_Com/N_Dist
## Slc39a5 2.070364 3.003691e-24 N_Prox/N_Com/N_Dist
## Ggt1 3.105118 3.576864e-24 N_Prox/N_Com/N_Dist
## Ldhb.7 2.332231 6.241133e-24 N_Prox/N_Com/N_Dist
## Pcbd1.3 2.043439 2.726802e-23 N_Prox/N_Com/N_Dist
## Cpn1.1 2.623152 3.193002e-23 N_Prox/N_Com/N_Dist
## Phgdh.8 2.884691 1.248214e-22 N_Prox/N_Com/N_Dist
## Lgals1.12 -2.742854 1.309490e-21 N_Prox/N_Com/N_Dist
## 0610005C13Rik 2.225278 1.584856e-21 N_Prox/N_Com/N_Dist
## Clu.4 2.073325 4.494623e-21 N_Prox/N_Com/N_Dist
## Alas2.14 -2.586118 6.791025e-21 N_Prox/N_Com/N_Dist
## Glud1.3 2.132607 2.785794e-20 N_Prox/N_Com/N_Dist
## Hba-a2.13 -5.321656 6.168150e-20 N_Prox/N_Com/N_Dist
## Ttc36 4.087390 2.159100e-19 N_Prox/N_Com/N_Dist
## Tmsb4x.8 -3.405344 2.436396e-19 N_Prox/N_Com/N_Dist
## Hbb-bs.14 -5.159188 3.769772e-19 N_Prox/N_Com/N_Dist
## Dcn.12 -2.169356 6.936446e-19 N_Prox/N_Com/N_Dist
## Cyb5a.5 2.061208 1.193316e-18 N_Prox/N_Com/N_Dist
## Col1a2.13 -2.067640 4.111726e-18 N_Prox/N_Com/N_Dist
## Aldob 3.945110 5.924123e-18 N_Prox/N_Com/N_Dist
## Slc34a1 3.340307 6.437418e-18 N_Prox/N_Com/N_Dist
## Hbb-bt.14 -5.166921 1.171618e-17 N_Prox/N_Com/N_Dist
## Hba-a1.13 -5.129147 3.504662e-17 N_Prox/N_Com/N_Dist
## Vim.10 -2.263959 2.836193e-16 N_Prox/N_Com/N_Dist
## Keg1 2.088222 2.979535e-16 N_Prox/N_Com/N_Dist
## Hmgcs1.2 2.009830 1.987327e-15 N_Prox/N_Com/N_Dist
## Spink1 2.402158 9.335577e-15 N_Prox/N_Com/N_Dist
## Sult1d1 3.400172 9.478032e-15 N_Prox/N_Com/N_Dist
## Mep1a 2.132565 1.937744e-14 N_Prox/N_Com/N_Dist
## Col3a1.13 -2.021109 2.335123e-14 N_Prox/N_Com/N_Dist
## Mest.12 -2.010056 3.614430e-14 N_Prox/N_Com/N_Dist
## Mkrn1.14 -2.139525 4.679515e-13 N_Prox/N_Com/N_Dist
## Tmem174 2.075951 2.329875e-12 N_Prox/N_Com/N_Dist
## Ass1.1 2.253869 4.666053e-12 N_Prox/N_Com/N_Dist
## Col1a1.13 -2.133876 5.454298e-12 N_Prox/N_Com/N_Dist
## Mt1.10 3.561380 2.923428e-10 N_Prox/N_Com/N_Dist
## Sostdc1.2 2.013489 9.531778e-10 N_Prox/N_Com/N_Dist
## Fabp3 2.024052 6.916508e-09 N_Prox/N_Com/N_Dist
## Cdkn1c.11 -2.265757 6.929272e-09 N_Prox/N_Com/N_Dist
## Hdc.1 2.094512 1.997294e-08 N_Prox/N_Com/N_Dist
## Mt2 2.434612 6.027492e-08 N_Prox/N_Com/N_Dist
## Spp2 2.217776 9.397472e-08 N_Prox/N_Com/N_Dist
## Hbb-y.12 -5.097419 1.151633e-04 N_Prox/N_Com/N_Dist
\(~\)
\(~\)
\(~\)
\(~\)